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The Heisenberg uncertainty principle for material objects 
is an essential corner stone of quantum mechanics and clearly 
visualizes the wave nature of matter. Here we report a 
demonstration of the Heisenberg uncertainty principle for the 
most massive, complex and hottest single object so far, the 
fullerene molecule C70 at a temperature of 900 K. We find a 
good quantitative agreement with the theoretical expectation: 
Ax x Ap = h, where Ax is the width of the restricting slit, 
Ap is the momentum transfer required to deflect the fullerene 
to the first interference minimum and h is Planck's quantum 
of action. 



Complementarity is one of the essential paradigms of 
quantum mechanics jD. Two quantities are mutually 
complementary in that complete (or partial) knowledge 
of one implies the complete (or partial) uncertainty about 
the other and vice versa. The most generally known case 
is the complementarity between position and momentum, 
as expressed quantitatively in the Heisenberg uncertainty 
principle Air x Ap > h/2. For neutrons the uncertainty 
relation has been demonstrated already back in 1966 by 
Shull Following the growing experimental efforts in 
atom optics during the last decade, the uncertainty prin- 
ciple has shown up implicitly in several experiments and 
has also been explicitly investigated in both the spatial 
and in the time domain g. 

While being a physical phenomenon of interest in its 
own right, the complementarity between momentum and 
position is also an important factor for practical pur- 
poses: for example it is applied for the preparation of 
transverse coherence in all experiments using collimated 
beams, a fact that can be mathematically phrased using 
the van Cittert-Zernike theorem |@-@|. 

There are good reasons to believe that complemen- 
tarity and the uncertainty relation will hold for all suffi- 
ciently well isolated objects of the physical world and that 
these quantum properties are generally only hidden by 
technical noise for larger objects. It is therefore interest- 
ing to see how far this quantum mechanical phenomenon 
can be experimentally extended to the macroscopic do- 
main. 

Here we report on an experiment investigating for the 
first time in a quantitative way the uncertainty relation 
upon diffraction at a single slit for a molecule as complex, 
massive and hot as the fullerene C70 (m — 840 amu) at 
an internal and translational temperature of 900 K. 

It is well known that the limit h/2 of the uncertainty 



relation Ax x Ap > h/2 is only reached for particular 
wave packets, for example of the Gaussian type. Evi- 
dently, the wave packet after passage through a rectan- 
gular slit is very different from this minimal uncertainty 
shape. This is also reflected in the far field distribution 
which is described by the well known sinc-function rather 
than a Gaussian. It is therefore a matter of definition 
and convenience which quantities to take as a measure 
of the position and momentum uncertainty in our case. 
Obviously, for a wave traversing a slit, one can take the 
slit width to be the measure of the spatial uncertainty 
Ax. The momentum uncertainty Ap can be related to 
the angular spread due to diffraction at the slit. Quanti- 
tatively, we define it as the momentum required to reach 
the first order interference minimum, which for small 
angles, as used in our experiment, lies in the direction 
6 ~ A/Ax = Ap/p = XAp/h. With this definition of the 
uncertainty we expect the following relation 

Ax x Ap = h 

to hold for the experiment. It is this equality that we 
quantitatively test in the following for a beam of C70. 
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FIG. 1. Setup of the experiment. A thermal C70 beam is 
produced by sublimation of fullerene powder at 900 K. The 
beam is narrowed by Si and diffracted by S2. Si is fixed to 
10 jj,m. The width of slit S2 is varied with ±30 nm accuracy 
for Ax < 1 fim. 



In the actual experiment the minima are not very well 
defined due to the broad de Broglie wave spectrum of 
our fullerene source. However, we can determine the full 
width at half maximum (FWHM) values of the central 
lobe with good precision. The far field diffraction pattern 
after a single slit is given by 1(9) ~ (sin(u)/u) 2 , where 
u = Ax • k6/2 and k — 2tt/X is the wave vector of the 
incident fullerene molecule. The half width at half maxi- 
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mum (HWHM) for the diffraction curve is then given by 
the solution of (sin(u 1 / 2 )/wi/2) 2 = 1/2, which yields 
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Practically, this means that we obtain the momentum 
uncertainty as defined above from the measured FWHM 
via Ap = Api/a/0.89. 

The setup of the experiment, shown in fig.l, is similar 
to that described in a previous publication Q. An effu- 
sive thermal fullerene beam is produced at about 900 K. 
The velocity spread was as large as Av/v ~ 0.6 and was 
taken into account into the numerical description of the 
experiment. 

The molecular beam is collimated by two piezo- 
controlled slits. The width of the first slit S%, is fixed 
at fO /jm, while the width Ax of the second slit, S-2 - 
which is located at the distance L\ — 113 cm further 
downstream - can be varied to investigate the position- 
momentum uncertainty relation. 

In order to also quantitatively describe the experiment 
the properties of the slits have to be known rather pre- 
cisely. The slits (Piezosysteme Jena) are made of two sili- 
con edges mounted on piezo-controlled flexure stages. We 
obtain information about the slit opening in three differ- 
ent ways: from the applied piezo voltage, from the read- 
ing of a strain gauge mounted to the slits and finally from 
the total number of molecules passing through the slit at 
a given opening. While the piezo voltage can be kept 
stable to better than AU/U < 10 -4 it is well known that 
piezos show creep, hysteresis and non-linearities. How- 
ever, it turned out in the experiments that the passive 
stability over a typical time of 1 hour was of the order of 
50 nm, as can be judged from the stability of the diffrac- 
tion patterns. From a calibration of the hysteresis curve 
we determine the change of the slit opening as a function 
of the piezo voltage. In order to know the absolute slit 
width we determined the zero position by measuring the 
number of molecules passing through the slit, when it was 
being closed. We estimate this method to be accurate to 
within ±30 nm. 

We extract the momentum spread Ap after S2 from 
the FWHM of the detected molecular beam AX exp in 
the detection plane, which is separated from S2 by the 
length L 2 = 133 cm. The observed distribution function 
fexp(x) = D(x) (8> M(x) is actually a convolution of the 
detector resolution function D{x) and the real molecular 
beam profile M{x). 

The scanning laser ionization detector has been char- 
acterized in depth in a previous publication |10|. For our 
present experiments with C70 the FWHM of the detector 
response was determined tobeZ? = 10±0.5 fim at a laser 
power of P=10 W. The effective FWHM detector height 
at this power was measured to be ~ 1 mm. 
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FIG. 2. Measured molecule distribution in the detection 
plane after passing through a piezo-controlled silicon slit hav- 
ing a width of Ax = 1.4 pan (top) and Ax ~ 70nm (bottom). 
Both the quantum mechanical calculation (continuous line) 
and the experiment (circle) show an increase of the beam 
width when going from medium (top) to narrow (bottom) slit 
widths. The dotted line indicates the wave calculation before 
the convolution with the known detector profile. 

The second contribution, related to the measured 
molecular beam profile M(x), is composed of both the 
classical collimation and the momentum spread due to 
the quantum uncertainty. In order to compare the ex- 
periment with the uncertainty relation derived before, 
we concentrate in the following on the half-width val- 
ues of these components only. Since the classical FWHM 
shadow width AX c i and the quantum contribution AX qu 
are completely independent their influence can be added 
quadratically to yield the FWHM value of M(x), which 
we denominate as AX mo i. The classical contribution 
AX c i can be derived from a simple geometrical shadow 
model. Taking the measured and the classically expected 
widths we can then deduce the contribution to the beam 
width due to the quantum uncertainty and we finally 
relate this spatial information to the corresponding mo- 
mentum uncertainty, which then reads: 



Ap = 
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where p z is the longitudinal momentum of the molecule. 

To trace out the uncertainty relation we varied the 
width of the second slit from about 20 fim down to 
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roughly 50 nm and record the molecular beam width 
in the detection plane. In figure 2 we show the mea- 
sured molecular beam profiles as full circles for two dif- 
ferent widths of the second collimation slit. We see a 
relatively narrow beam of AX exp = 17 /xm for the slit 
width Ax = 1.4 /im (fig 2a) and again a strong growth to 
AX exp = 43 fim for the slit width Ax = (0.07 ±0.03 )fj,m 
(fig 2b). The error bars in fig. 2 represent the statistical 
uncertainty due to the very low count rate, in particular 
at the smallest slit width. 

The dashed line follows a full wave calculation as de- 
scribed below in order to show the molecular beam profile 
as given by diffraction alone. The continuous curves rep- 
resent the same model but convoluted with the detector 
profile. 

Close inspection of the data shows a good agreement 
between the convoluted wave model and the experimental 
data. This good agreement is the first demonstration of 
single slit diffraction for a molecule as heavy, complex 
and hot as C70. 

Fig. 2b is actually an interesting complement to 
the high contrast interference fringes of fullerenes after 
diffraction at a nanofabricated grating with a grating 
constant of 100 nm, which we could demonstrate in a 
previous publication 0|. The single slit pattern shown 
here is the envelope of the far-field grating interference 
pattern. This provides a striking proof of the wave na- 
ture of the fullerene C70 because it demonstrates that 
the previous minima must have been due to destructive 
interference. 
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FIG. 3. The experimental molecular beam width (full cir- 
cle) is compared with the quantum prediction (continuous 
line) as a function of the slit opening Ax. The agreement is ex- 
cellent across the whole range of slit openings (70nm...20 /im). 
A purely classical shadow model predicts the dotted line and 
is in marked disagreement with the data for Ax < 4/xm. The 
latter is therefore designated as the quantum regime and mag- 
nified in the inset of fig 3. 

From the whole series of experiments with varying slit 
widths we have extracted the FWHM values from the 
experiment and we compare them with a quantum wave 



model in fig. 3. An excellent agreement between ex- 
pectation and experiment is found throughout the whole 
range of values. We can distinguish essentially two dif- 
ferent regimes corresponding to a pure quantum regime 
(left part of fig. 3) and a range which can be very well 
described using a classical ball model (right part of fig. 
3). The continuous wave calculation curve and the dot- 
ted classical line coincide almost completely down to a 
slit width of about Ax = 4 /im. Below this value the 
quantum mechanical momentum spread Ap contributes 
significantly to the beam width in the detection plane. 
This quantum range is magnified in the inset of fig. 3. 

The horizontal error bars in this picture have two com- 
ponents, namely the precision both of the absolute zero 
and of the scaling of the piezo translation as a function of 
the applied voltage. Both are only important for a small 
slit width Ax. The absolute zero (closed) position of the 
piezo slits is known with an error of ± 30 nm, as men- 
tioned above. The scaling with the applied piezo voltage 
is non-linear and follows a hysteresis curve, which has 
been calibrated. We estimate an uncertainty of ±3% in 
the calibration of the hysteresis curve. 

The vertical error bars estimate the uncertainty of the 
measured width of the beam in the detection plane. For 
small slit widths Ax these values are obtained from a 
least squares Gaussian fit to the detected curve. For 
large Ax the marked trapezoidal shape as well as the 
high signal-to-noise ratio permit a direct reading of the 
experimental and theoretical FWHM values with very 
high accuracy. 

The numerical simulations in fig. 3 are based on the 
fact that the Schrodinger equation of our time indepen- 
dent problem is formally equivalent to the Hclmholtz 
equation and can therefore be treated using all the meth- 
ods well known from optics. The solution is done in close 
analogy to the numerical approach as used in [jl2| for 
neutrons and similarly in [pL3|Jl4| l for atoms. There is no 
free parameter in the calculation except for a broadening 
of the detector resolution by 3.5/xra with respect to the 
best detector resolution curves recorded some time ear- 
lier (nj . This offset is most likely explained by a residual 
tilt of 2.7 mrad between laser and diffraction slit. This 
is in agreement with diffraction curves not shown here, 
which were recorded using the same setup but at half the 
width of the first collimation slit. 

Since in previous papers it has been pointed out that 
the form factor of single slit diffraction may be influenced 
by the van der Waals interaction between the molecule 
and the slit walls ^|) one may wonder whether this ef- 
fect may become visible in the present experiment. How- 
ever, the slit widths here, except for the smallest, are 
much larger than in the former grating diffraction ex- 
periments, where the effective slit width was reduced by 
about 15 nm. Since the van der Waals potential above a 
surface decreases with the third power of the object- wall 
distance the effect becomes small for the present study 
although the slit thickness is bigger than that of the pre- 
viously used SiN x gratings. For the smallest slit width, 
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Ax ~ 70nm, a possible contribution is masked by the 
experimental error bar. 

In the following we compare our findings with the 
Heisenberg uncertainty relation between position and 
momentum. For this we use the method as indicated fur- 
ther above: From the measured beam width we separate 
the influence of the detector resolution in a deconvolu- 
tion procedure. The remaining molecular beam width is 
then decomposed into its classical and quantum part. 
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FIG. 4. Experimental verification of the uncertainty rela- 
tion for C70. The momentum values are derived from the far 
field molecular beam widths as described in the text. The po- 
sition uncertainty is given by the width of the second slit. The 
experiment is in very good agreement with the wave model 
expectation: Ax x Ap = h and in marked disagreement with 
a purely classical shadow model. 

We can then plot Ap as derived from eqn. 1 as a 
function of Ax for slit openings lying well in the quan- 
tum regime and obtain fig. 4. The full circles repre- 
sent the values extracted from the experiment with error 
bars directly related to those of the inset of fig. 3. The 
continuous line follows our original expectation function 
Ax x Ap — h without any additional fit parameter. We 
regard the good quantitative agreement between the data 
points and the predicted curve as a good support for the 
validity of the Heisenberg uncertainty principle for the 
fullerene C70 - i.e. in a complexity region not studied so 
far. 
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